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Abstract 

One measure of geometrical complexity of a spatial curve is the number of crossings in a planar 
projection of the curve. For iV-noded ring polymers with a fixed knot type, we evaluate numerically 
the average of the crossing number over some directions. We find that the average crossing number 
under the topological constraint are less than that of no topological constraint for large N. The 
decrease of the geometrical complexity is significant when the thickness of polymers is small. The 
simulation with or without a topological constraint also shows that the average crossing number 
and the average size of ring polymers are independent measures of conformational complexity. 



PACS numbers: 05.40.Fb, 36.20.-r, 61.41. +e 



I. INTRODUCTION 



Complexity of conformations of polymer chains should play an important role in the 
physics of polymers. Various dynamical or statistical properties of entangled polymers in 
solutions or melts should be related to the complexity of polymer chains as space curves. For 
instance, the reptation theory takes into account the entanglement complexity of polymer 
chains around a polymer by assuming the tube along it 

However, it is not trivial to investigate any aspect of the complexity of conformations of 
polymers directly through computer simulation. In fact, it is not clear how to express the 
entanglement complexity numerically for mutually entangled polymers. Furthermore, it is 
not known how the complexity of polymer conformations should depend on the property 
that polymer chains can not cross each other. 

In this paper, we discuss the geometrical complexity of conformations of a ring polymer 
under some topological constraint. As a measure of geometrical complexity of a spatial 
curve, we consider the number of crossings in a planar projection of the curve, and take the 
average over some directions 0, |3|. Through computer simluation, we evaluate the average 
crossing number for iV-noded ring poymers with fixed knots, and discuss their behaviors 
with respect to N. We thus make it clear how the topological constraints can modify the 
geometrical complexity. Here we note that the condition that a ring polymer never cross 
with itself corresponds to a topological constraint on it, as far as its statistical properties are 
concerned. Furthermore, we plot the graphs of the mean square radius of gyration versus 
the average crossing number, for the SAPs with fixed knots. Then we find that the average 
crossing number and the mean-squared gyration radius behave independently under the 
topological constraints. It is thus suggested that the relation between the average sizes of 
ring polymers and the geometrical complexity of the conformations should be non-trivial. 

Let us explain the average crossing number, more precisely. The writhe of a linear or 
ring polymer is defined by the average of the number of signed crossings appearing in a 
projection of the curve over all directions. As a simplified version of the writhe, Janse 
van Rensburg et.al introduced the number of crossings 0. It is defined by the number of 
unsigned crossings in a projection. Then, the average number of crossings is defined by the 
number of crossings averaged over all possible projections. There are several numerical or 
theoretical studies on the average crossing number [|, @, [jj . 
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The average crossing number is also related to ideal knots. For a knot, the ideal knot 
is given by its tightest geometric configuration P, PH|. Katritch et. al have obtained 
the average crossing number for the ideal knots up to 11 essential crossing numbers. The 
average crossing number should be useful for flexible DNA knots in thermal equilibrium ||. 
Furthermore, it should also be useful for statistical or dynamical studies on knotted ring 
polymers @, 0, . 



II. METHODS OF SIMULATION 



We consider self-avoiding polygons (SAPs) consisting of N rigid impenetrable cylinders 
of unit length and radius r. There is no overlap allowed between any pair of non-adjacent 
cylindrical segments, while next-neighboring cylinders may overlap each other. We call them 
cylindrical SAPs, for short. A large number of cylindrical SAPs can be constructed by the 



cylindrical ring-dimerization method [|13|]. The method is based on the algorithm of ring 



dimerization |TJ|], and useful for generating long SAPs systematically. In this paper, we 
have constructed M=10 4 samples of N cylindrical segments with radius r, where N is from 
20 to 1000 and r is from 0.0 to 0.07. We determine the number Mk of polygons with a knot 
K, enumerating such polygons that have the same set of values of the two invariants: the 
determinant of the knot A/<(— 1) and the Vassiliev invariant V2(K) of the second degree. 

The mean value A of the average crossing number for a set of SAPs is defined by A = 
'Yl i=1 Ai/M . Here, A4 denotes the average crossing number of the ith polygon. The mean 
value of the average crossing number for a set of the SAPs with a knot K is given by 
Ak = J2i=i Ak,i/Mk, where Ak,% denotes the average crossing number of the ith SAP 
having the knot K. Thus, A is given by the average of Ax 's over all knots. 

In the simulation, we have obtained the average crossing number A and A#-'s for the 
trivial and trefoil knots in the range of N from 20 to 1000. Here, we evaluate the average 
crossing number by taking the average over the x, y and z directions. 

The mean square radius of gyration R 2 for SAPs with N nodes is given by R 2 = 
Y^n '771=1 (C^n — -Rm) 2 )/2A 2 . Here R n is the position vector of the nth node, and the symbol 
(■) denotes the average over M polygons generated. For a knot K, we define the mean square 
radius of gyration R\ for SAPs of the knot K by R 2 K = R 2 Ki jM^ p, 0. Here, R 2 Ki 

denotes the mean square radius of gyration for the i th SAP in the Mk polygons of the knot 
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K. In terms of -R^'s, R 2 is given by the average over all knots: R 2 = J2k M k R 2 k /M. 

III. AVERAGE CROSSING NUMBER UNDER A TOPOLOGICAL CON- 
STRAINT 

We now discuss the N dependence of the average crossing number under the topo- 
logical constraint of a knot K. The double-logarithmic graph of A K versus iV is given in 
Fig. 1 for the cylindrical SAPs of cylinder radius r = 0.003, for the trivial and trefoil knots. 
The graphs can be approximated by some straight lines in some large N region. 

Let us consider the ratio of A K to A for the two knots. In Fig. 2, the ratio A K /A versus 
N is plotted in a double logarithmic scale for the trivial and trefoil knots. The graph of 
the trivial knot has a concave curve: the ratio A triv /A is almost constant with respect to N 
for small N and then decreases with a larger gradient for large N. On the other hand, the 
graph of the trefoil knot suggests that A tre j A could be roughly approximated by a power of 
N for some finite values of N. 

In Fig. 2, the ratios A K /A's for the two knots become less than 1.0 when N is large, i.e., 
the Ak$ are less than the average crossing number A averaged over all knots for large N. 
Thus, the topological constraints make the conformations of ring polymers simpler with re- 
spect to the geometric complexity, when N is large enough. The reduction of conformational 
complexity may be related to entropic repulsion arising from the topological constraints. 

For the trivial knot, A triv /A is smaller than 1.0 in the whole range of N, while for the 
trefoil knot A tre /A is larger than 1.0 for N < 200, and smaller than 1.0 for N > 300. Thus, 
with respect to the average crossing number, the iV-noded ring polymers with the trivial 
knot are less complex for any N than those of no topological constraint , while those of the 
trefoil knot are more complex for small N and less complex for large N than those of no 
topological constraint. 

The average crossing number of the trivial knot, A triv , is smaller than that of the trefoil 
knot, A tre , through the whole range of N. It suggests that the more complicated knot should 
have the larger value of the average crossing number, at least for finite N. The tendency 
is also seen in the data for the different values of r. As N increases, however, the average 
crossing numbers A^'s of the two knots gradually become close to each other and they 
become almost the same value when N is very large. 
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The average crossing number Ak of a knot K have similar properties with the inverse of 
the mean-squared gyration radius R 2 K of ring polymers with the knot K. In our previous 
work [|15|, |I6[| , it is shown that for some random polygons and cylindrical SAPs, the double- 
logarithmic graph of the ratio R 2 K / R 2 versus N is given by a downward convex curve and is 
larger than 1.0 for the trivial knot, while for the trefoil knot the graph is given by a straight 
line and the ratio R 2 K /R 2 is smaller than 1.0 for small N and larger than 1.0 for large N. 

We now consider possible asymptotic behaviors of the average crossing number. Let us 
review some known results on the large N behavior of iV-step self-avoiding walks (SAWs). 
The average crossing number can also be defined for SAWs, and we denote it by A$aw ■ 
In Ref. 0, it is discussed that for asymptotically large N, the mean value A S aw of the 
average crossing number of self-avoiding walks (SAWs) is given by Agj^/N = a — bN 1 ^ 2 " 
with some constants a and b. Here v is given by the exponent of the average size of SAWs. 
We also note that in Ref. ||, the large N behavior of A$aw is approximated by a power of 
N, A SAW ~ N^ SAW , with the effective exponent: /i SAW =1-122 ± 0.005. 

In order to illustrate some r dependent properties of Ak, let us introduce an asymptoic 
expansion for the ratio Ak/A versus N. Based on the asymptotic expansion of Asaw i n 
Ref. 0, we assume the asymptotic expansion with repsect to iV as follows: Ak/N = cik — 
bK N 1 ~ 2uk for any knot K. Here and 6^ are fitting parameters. From the simulations fT5] , 
|To| , we may assume that vk = v- For the ratio Ak/A, we thus have the following formula: 
Ak/A = oik (1 — Pk N^ 2u ). Here, and (3k are fitting parameters corresponding to a K /a 
and bK/ctK — b/a, respectively. We apply it to the data of Ak/A of the ./V-noded SAPs with 
TV larger than some cut-off and for the different values of r. 

Let us discuss the best estimates of a K /a plotted against cylinder radius r in Fig. 3. 
When r is small, the ratio a^/a becomes smaller than 1.0 both for the trivial and trefoil 
knots. This is consistent with the observation of Fig. 2 that the ratio Ak/A for r = 0.003 
decreases against iV and is smaller than 1.0 for large N. In Fig. 3, the ratio a^/a increases 
monotonically with respect to cylindrical radius r, and it becomes close to the value 1.0 at 
some large value of r . If the ratio Ak/A becomes 1.0, then there is no topological effect 
on the average crossing number. On the other hand, if Ak/A is less than 1.0, then it may 
be a consequence of the topological constraint. Thus, the behavior that the graph of a^/a 
versus r increases upto 1.0 suggests that topological constraints on ring polymers make their 
conformations simpler for small r, while the topological effect becomes weak for large r. 
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Let us consider again the r dependence of Ak that the graph of a^/a versus r increases 
upto 1.0 for the two knots, as shown in Fig. 3. A very similar behavior has also been observed 
for the case of the ratio R 2 K /R 2 of the mean squared gyration radii of the cylindrical SAPs 
with radius r. In Refs. ]T5| , it has been shown for the cylindrical SAPs that a topological 
constraint on a ring polymer gives effective expansion to it, i.e., the ratio R 2 K /R 2 becomes 
larger than 1.0 for large N. Furthermore, the effective expansion is significant when the 
cylinder radius r is small; the large N limit of R\jR 2 decreases to 1.0 as r increases. 

Summarizing the simulation results, we may conclude that the effect of topological con- 
straints should be significant when the ring polymer is thin, both for the average crossing 
number A K and the mean square radius of gyration R\. 

IV. THE MEAN SQUARED GYRATION RADIUS AND THE GEOMETRIC 
COMPLEXITY 

Let us discuss that the relation between the average crossing number and the average 
size of ring polymers should be nontrivial. In Fig. 4, the mean square radii of gyration 
R 2 and Rj^'s are plotted against the average crossing numbers A and i^'s in a double 
logarithmic scale for the cylindrical SAPs of r=0.003. The estimates of R 2 and -R|-'s in Fig. 
4 are different for any given value of the average crossing numbers. Thus, it follows that the 
mean-squared gyration radius and the average crossing number are independent quantities 
describing some geometric properties of conformations of ring polymers. 

With the same average crossing number given, R 2 riv and R 2 re are larger than R 2 for the 
cylindrical SAPs, as shown in Fig. 4. Thus, we may say that the topological constraints 
make the average sizes of ring polymers larger with respect to the average crossing number. 
Among the -R^'s, R\ riv is larger than R 2 re in Fig. 4. It is thus suggested that the more 
complicated knot should have the smaller radius of gyration, as far as some finite values of 
Ak's are concerned. 

When the average crossing number is very large, the mean-squared gyration radius R 2 
(or R 2 K ) can be approximated by some power of A (or Ak): R 2 ~ A 1 for the average over all 
knots; R 2 K ~ AJf for the trivial and trefoil knots. Applying the power law approximation 
to the data of N > 100, we obtain the following estimates of the effective exponents: 7 = 
0.860 ± 0.001, 7 trfo = 0.952 ± 0.004 and ^ tre = 1.069 ± 0.006. The graphs of Fig 4 suggests 
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that the effective exponent 7^ should be independent of the knot type. Similar results are 
obtained also for the SAPs with the different values of r. 



Acknowledgments 

We would like to thank Prof. K. Ito for helpful discussions. 



M. Doi and S.F. Edwards, The Theory of Polymer Dynamics (Clarendon Press, Oxford, 1986). 
E. J. Janse van Rensburg, D. A. Sumners, E. Wasserman and S. G. Whittington, J. Phys. A: 
Math. Gen. 25 (1992) 6557. 

E.Orlandini et.al, J.phys.A: Math Gen. 27 (1994)L333. 
G.A. Arteca, Phys. Rev. E 49 (1994) 2417. 
J-Y Huang and P-Y Lai, Phys. Rev E 63 (2001) 021506. 
A.L. Kholodenko and T.A. Vilgis, Phys. Rep. 298 (1998) 251. 
P. Grassberger, J. Phys. A : Math. Gen. 34 (2001) 9959. 

V. Katritch, J. Bednar, D. Michoud, R. G. Scharein, J. Dubochet, and A. Stasiak, Nature 
(London) 384, 142 (1996) 

A. Stasiak, J. Dubochet, V. Katritch, and P. Pieranski, in Ideal Knots edited by A. Stasiak, 

V. Katritch and L. H. Kauffman (World Scientific, Singapore, 1998) pp. 1-19. 

A. Yu. Grosberg, A. Feigel, and Y. Rabin, Phys. Rev. E 54, 6618 (1996). 

A. Yu. Grosberg, in Ideal Knots edited by A. Stasiak, V. Katritch and L. H. Kauffman (World 

Scientific, Singapore, 1998) pp. 129-142. 

M.K. Shimamura and T. Deguchi, Phys. Rev. E 66, R040801(2002). 
M.K. Shimamura and T. Deguchi, Phys. Lett. A 274, 184 (2000). 

Y.D. Chen, J. Chem. Phys. 74, 2034 (1981) ; J. Chem. Phys. 75, 2447 (1981); J. Chem. Phys. 
75, 5160 (1981). 

M.K. Shimamura and T. Deguchi, Phys. Rev. E 64, R 020801 (2001). 
M.K. Shimamura and T. Deguchi, Phys. Rev. E 65, 051802 (2002). 



7 



FIG. 1: The average number of crossings with a knot type K Ax versus N for r=0.003. Nu- 
merical estimates of Ak for i^=trivial and trefil knots are shown by closed triangles and squares, 
respectively. 

FIG. 2: Double- logarithmic plots of the ratio Ak/A versus N for cylindrical SAPs with r=0.003. 
A triv /A and A tre /A are shown by closed triangles and squares, respectively. 

FIG. 4: The mean-square radius of gyration R K versus the average number of crossings Ak of 
cylindrical SAPs for r=0.003. Numerical estimates for i"C=trivial knot, 3i are shown by closed 
triangles and squares, respectively. Data of R 2 versus A are shown by closed circles. N are given 
by 51, 151 and 100j+l for j = l,- • • ,10. 



FIG. 3: The ratio clk /a versus cylinder radius r for cylindrical SAPs. The values of atriv/a and 
atre/a are shown by closed triangles and squares, respectively. 
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